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We study the heavy-quark contributions to the proton structure functions F2(x,Q^) and 
F2(x,Q^), with i = c, 6, for small values of Bjorken's x variable and provide compact 
formulas for their ratios Ri = FI/F2 that are useful to extract F2(x,Q^) from measure- 
ments of the doubly differential cross section of inclusive deep- inelastic scattering at DESY 
HERA. Our approach naturally explains why Ri is approximately independent of x and 
the details of the parton distribution functions in the low-x regime. 

1 Introduction 

The totally inclusive cross section of deep-inelastic lepton-proton scattering (DIS) depends on 
the square s of the centre-of-mass energy, Bjorken's variable x = Q'^/{2pq), and the inelasticity 
variable y = / {xs), where p and q are the four- momenta of the proton and the virtual photon, 
respectively, and Q'^ = —q'^ > 0. The doubly differential cross section is parameterized in terms 
of the structure function F2 and the longitudinal structure function Fl^ as 

^^ = '-^ni + ii-,m(^.Q')-i-F.M')}. (I) 

where a is Sommerfeld's fine-structure constant. At small values of x, Fl becomes non-negligible 
and its contribution should be properly taken into account when the F2 is extracted from the 
measured cross section. The same is true also for the contributions F2 and F£ of F2 and Fl 
due to the heavy quarks i = c^b. 

Recently, the HI [H [2 [3] and ZEUS H El [6] Collaborations at HERA presented new data 
on F2 and F2. At small x values, of order 10~^, F2 was found to be around 25% of which 
is considerably larger than what was observed by the European Muon Collaboration (EMC) at 
CERN [7J at larger x values, where it was only around 1% of F2. Extensive theoretical analyses 
in recent years have generally served to establish that the data can be described through 
the perturbative generation of charm within QCD (see, for example, the review in Ref. [8j and 
references cited therein). 

In the framework of Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) dynamics [9], 
there are two basic methods to study heavy- flavour physics. One of them [TO] is based on 
the massless evolution of parton distribution functions (PDF) and the other one on the photon- 
gluon fusion (PGF) process [11]. There are also some interpolating schemes (see Ref. [12] and 
references cited therein). The present HERA data on [2l[3l[5l[6] are in good agreement with 
the modern theoretical predictions. 
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In earlier HERA analyses [H S] , F£ and were taken to be zero for simplicity. Four years 
ago, the situation changed: in the papers [21 [3l [3 [6] , the F£ contribution at next-to- leading 
order (NLO) was subtracted from the data. 

In this paper, we present compact \ow-x approximation formulae [l3j for the ratio Ri = 
FI/F2 at leading order (LO) and NLO, which greatly simplify the extraction of from mea- 
surements of dp a^'^ / {dx dy) . 



2 Parton distribution functions at small x 

The standard program to study the small x behavior of quarks and gluons is carried out 
by comparison of the data with the numerical solution of the DGLAP equations fitting the 
parameters of the x profile of partons at some initial Qq QCD energy scale A (see, for 

instance. Oils]). However, in analyzing exclusively the small x region {x < 0.01), there is 
the alternative of doing a simpler analysis by using some of the existing analytical solutions 
of DGLAP in the small x limit (see [16] for review). It was done in Refs. [I7]-[IS]5 where 
it was pointed out that the HERA small x data can be interpreted in the so called doubled 
asymptotic scaling (DAS) approximation related to the asymptotic behavior of the DGLAP 
evolution discovered in [20j many years ago. 

Here we illustrate results obtained in [TSl [19]: the small x asymptotic PDF form in the 
framework of the DGLAP equation starting at some Qq with the flat function: 

xfa{x,Ql) = Aa (hereafter a = g,^), (2) 



where xfa are the leading-twist PDF parts and Aa are unknown parameters that have to be 
determined from data. We neglect the non-singlet quark component at small x. 

We would like to note that HERA data [21] show a rise of F2 at low values (Q^ < IGeV^) 
when X ^ 0. This rise can be explained naturally by incorporation of higher-twist terms in the 
analysis (see [19] and Fig.l). 

We shortly compile the LO results (the NLO ones may be found in [IHIIIS]), which are: 

/a(x,Q2) = f+{x,Q^) + f-{x,Q^), (3) 
ftM^) = (A, + ^A^h{a)e-~'+^'> + 0{p) , (4) 

f^{x,Q') = {(^A, + ^A,)p7i(a)e-^+«^ + 0(p), (5) 

/-(x,Q2) = + o{x), (6) 

f-{x,Q^) = A^e-'^-W^ + 0{x), (7) 

where where e = (Xli ^f)// is the average charge square and (i+(l) = 1 + 20//(27/3o) and 
d-{l) = 16//(27Ad) are the regular parts of d-^ and d- anomalous dimensions, respectively, in 
the limit n ^ 1 lj. The functions 1^ = 0^1) are the modified Bessel functions 1^ and the 

^ For a quantity k{n) we use the notation k{n) for the singular part when n ^ 1 and k{n) for the corresponding 
regular part. 
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Figure 1: The structure function F2 as a function of x for different bins. The sohd and 
dashed hues are obtained without and with higher- twist terms, respectively. 



variables a and p are given by 

A 1 / N / d-i-S <J ; 12 

<x = 2Vd+.ln(x), p=y^ = -^-^, d+ = -^, (8) 

where Po is the first coefficient of the QCD beta function and s = ln[as(Qo)/<^s(Q)], with Qo 
being the initial scale of the DGLAP evolution, and as{p) = as{fi)/{4:7r) is the couplant with 
the renormalization scale /i. 



3 Master formula 

We now derive our master formula for Ri{x,Q'^) appropriate for small values of x, which has 
the advantage of being independent of the PDFs fa{x^ Q^). In the low-x range, where only the 
gluon and quark-singlet contributions matter, while the non-singlet contributions are negligibly 
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small, we hav^ 

F^(x,Q2)= ^ Cijx,Q')^xfg{x,Q'), (9) 
^=+,- 

where / = ± labels the usual + and — linear combinations of the gluon contributions, Cl ^(x, Q^) 
are the DIS coefficient functions, which can be calculated perturbatively in the parton model of 
QCD, and the symbol denotes convolution according to the usual prescription, f{x) ®g{x) = 
{dy / y) f {y)g{x / y) . Massive kinematics requires that C^. ^ = for x > bi = 1/(1 + 4ai), where 
ai = ml/Q'^. We take to be the solution of mi{mi) = m^, where rni{fi) is defined in the 
modified minimal-subtraction (MS) scheme. 

Exploiting the low-x asymptotic behaviour of /^(x, Q^) [22] , 

1 
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.1+5, 



(10) 



where the rise of /^(x,(5^) as x ^ is less than any power of x, Eq. (|9]) can be rewritten as 
[231 [24] 

Flix^Q')^ Mijl^5i,Q')xfl{x,Q'), (11) 

where ^ 

MUn.Q') = dxx^-'Cijx^Q') (12) 
Jo 

is the Mellin transform, which is to be analytically continued from integer values n to real values 
1 + ^^ [25]. 

In the DAS approach!, one has M+^(1,Q2) = M-^(1,Q2) if M[^^{n,Q'^) are devoid of 
singularities in the limit 5i ^0, as we assume for the time being. Such singularities actually 
occur at NLO, leading to modifications to be discussed in Section [5l Defining M^^aO^iQ^) = 
M^^(l, Q^) and using (|9]), Eq. ([TT]) may be simplified to become 



F^(x,Q2)^M,,,(1,Q2)x/,(x,Q^). 
In fact, the non-perturbative input fg{x, Q^) does cancels in the ratio 



M2,<,(1,Q2)' 



(13) 



(14) 



which is very useful for practical applications. Through NLO, Mk,g{l, Q^) exhibits the structure 



Mk,g{l,Q'') = efas{^l)■{Mi°l{l,ai) + a,{^l) 



M«(l,a,) + M^;;(l,a,)ln^ 



(15) 



where Inserting Eq. (fT5|) into Eq. ([T4|) . we arrive at our master formula 



M^°l{l,a,)+a,{p) 




+ MgJ(l,ai)ln(MV"^f/ 


M^°^g{l,ai)+asiii) 


Mg(l,a,) 


+ Mg(l,a.)ln(/xVmf)" 



+ 0{al). (16) 



^Here and in the following, we suppress the variables /i and rrii in the argument lists of the structure and 
coefficient functions for the ease of notation. Moreover, a further simplification is obtained by neglecting the 
contributions due to incoming light quarks and antiquarks in Eq. (|9]), which is justified because they vanish at 
LO and are numerically suppressed at NLO for small values of x. One is thus left with the PGF contribution. 

^The singular PDF behavior has been considered recently in j26) . 
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We observe that the right-hand side of Eq. (p!6|) is approximately independent of a remark- 
able feature that is automatically exposed by our procedure. In the next two sections, we 
present compact analytic results for the LO (j = 0) and NLO (j = 1, 2) coefficients M^-^^(l, a^), 
respectively. 



4 LO results 

The LO coefficient functions of PGF can be obtained from the QED case [27] by adjusting 
coupling constants and colour factors, and they read [28l [29] 

a) = -2x{[l - 4x{2 - a)(l - x)]f3 - [1 - 2x{l - 2a) + 2x^(1 - 6a - 4a^)]L(/3)}, 
= Sx^[{l-x)(3-2axL{(3)], (17) 

where 



Performing the Mellin transformation in Eq. (fT2|) . we find (see details in [13j) 

= |[1 + 2(1- a) J(a)], Mi°) (1, a) = ^5[1 + 6a -4a(l + 3a) J(a)], (19) 

where 

J(a) = -v^lnt, t=^—^, (20) 

1 + v^ 

At LO, the low- a; approximation formula thus reads 

^ l + 6a.-4a.(l + 3a.)J(a.) 

1 + 2(1 - ai)J(ai) ^ ^ 

5 NLO results 

The NLO coefficient functions of PGF are rather lengthy and not published in print; they are 
only available as computer codes For the purpose of this letter, it is sufficient to work in 
the high-energy regime, defined by x <C 1, where they assume the compact form ^ 

Ci%,a)=/3Ri^l{l,a), (22) 

with 

4^](l,a) = ^CA[5 + (13-10a)J(a) + 6(l-a)/(a)], i?g(l, a) = -4CaM^;^](1, a), 

1(1. = -^CAb{l - 12a - [3 + 4a(l - 6a)]J{a) + 12a(l + 3a)/(a)}, (23) 
9 

where Ca = N for the colour gauge group SU(N), J (a) is defined by Eq. ([2Q|) . and 



I{a) = -Vb 



C(2) + i In^ t - ln(a6) Int + 2Li2(-t) 



(24) 
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Here, C(2) = 7r^/6 and Li2{x) = — fQ{dy/y) ln(l — xy) is the dilogarithmic function. 

As already mentioned in Section [3l the Mehin transforms of cj^g{x,a) exhibit singularities 
in the limit 6i 0, which lead to modifications in our formalism, namely in Eqs. (|13p and (fT6|) . 
As was shown in Refs. [231 [THl HHj, the terms involving 1/(5; depend on the exact form of the 
subasymptotic low-x behaviour encoded in fi{x, Q^), as 



1 

si 



fUy,Q'), 



where x = x/b. In the generahzed DAS regime, given by Eqs. ([3])-([7j), we have 



1 



1 hjajx)) ^ 
p Io{a{x))' 



0- X 



(25) 



(26) 



Because the ratio fg {x^Q'^)/ f^{x,Q'^) is rather small at the values considered, Eq. (p!3|) is 
modified to become 

FUx,Q^) « MkAl,Q^)xfg{x,Q^), (27) 

where Mfe^^(l,(5^) is obtained from Mk^gin^Q'^) by taking the limit n 
l/(n — 1) 1/^^+- Consequently, one needs to substitute 



Mg(l,a)^M^^^(l,a) (i = l,2) 
in the NLO part of Eq. (p!6|) . Using the identity 
1 



1 and replacing 
(28) 



U_^fx_ 

y \y 



h{<y{y)) 



h{a{x)) 

we find the Mellin transform ^ of Eq. ^ to be 



1 



ln(a6) 



J{a) 



1 



ln(a6) 



J{a) 



i?g^(l,a) (j = l,2). 



(29) 



(30) 



The rise of the NLO terms as x ^ is in agreement with earlier investigations 



6 Results 

As for our input parameters, we choose Qq = 0.306 GeV^, rric = 1.25 GeV and rrib = 
4.7 GeV. While the LO result for Ri in Eq. ([2T]) is independent of the unphysical mass scale 
/i, the NLO formula (p!6|) does depend on it, due to an incomplete compensation of the /i 
dependence of as{p) by the terms proportional to ln(/i^/(5^), the residual /i dependence being 
formally beyond NLO. In order to estimate the theoretical uncertainty resulting from this, in 
[13] we put /i^ = and vary ^. Besides our default choice ^ = 1 + 4a^, we also considered 
the extreme choice ^ = 100, which is motivated by the observation that NLO corrections are 
usually large and negative at small x values [33j . A large ^ value is also advocated in Ref. [34] , 
where the choice = 1/x^, with 0.5 < A < 1, is proposed. 

We now extract F2{x^Q'^) {i = c, 6) from the HI measurements of the cross sections in 
Eq. O at low (12 < < 60 GeV^) ^ and high (Q^ > 150 GeV^) ^ values of using 
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X 


F^{x,Q')-W (HI) 




F^{x 


,Q^).10^ (HI) 


F^{x 


g^)-io^ 


12 


0.197 


435 ± 78 


431 




45 ±27 




45 


12 


0.800 


186 ± 24 


185 




48 ±22 




48 


25 


0.500 


331 ± 43 


329 




123 ± 38 




123 


25 


2.000 


212 ±21 


212 




61 ±24 




61 


60 


2.000 


369 ± 40 


368 




190 ± 55 




190 


60 


5.000 


201 ± 24 


200 




130 ± 47 




130 


200 


0.500 


202 ± 46 


202 




413 ± 128 




400 


200 


1.300 


131 ±32 


130 




214 ± 79 




212 


650 


1.300 


213 ±57 


214 




243 ± 124 




238 


650 


3.200 


92 ±28 


91 




125 ±55 




125 



Table 1: Values of F2^(x, Q^) and F2^(x, Q^) extracted from the HI measurements of d^^ and d^^ 
at low [3] and high ^ values of (in GeV^) at various values of x (in units of 10~^) using our 
approach at NLO for ji^ = with ^ = 1 ± 4ac. The LO results agree with the NLO results 
within the accuracy of this table. For comparison, also the results determined in Refs. [21 [3] 
are quoted. 



the LO and NLO results for Ri derived in Sections 2] and [5l respectively. Our NLO results for 
/i^ = with = 1 ± Aai are presented for i = c, 6 in Table [TJ where they are compared with 
the values determined by HI. We refrain from showing our results for other popular choices, 
such as /i^ = 4mf,(5^ and even /i^ = lOOQ^ because they are very similar. We observe that 
the theoretical uncertainty related to the freedom in the choice of /i is negligibly small and find 
good agreement with the results obtained by the HI Collaboration using a more accurate, but 
rather cumbersome procedure [2l [3] . 

In order to assess the significance of and the theoretical uncertainty in the NLO corrections 
to we show in Fig. [2] the dependences of Rc, Rb^ and Rt evaluated at LO from Eq. ([2T]) 
and at NLO from Eq. ([16]) with /i^ = 4mf , ± 4mf . We observe from Fig. [2] that the NLO 
predictions are rather stable under scale variations and practically coincide with the LO ones 
in the lower regime. On the other hand, for 4mf , the NLO predictions overshoot the 

LO ones and exhibit a strong scale dependence. We encounter the notion that the fixed-flavour- 
number scheme used here for convenience is bound to break down in the large-Q^ regime due to 
unresummed large logarithms of the form ln((5^/m|). In our case, such logarithms do appear 
linearly at LO and quadratically at NLO. In the standard massless factorization, such terms are 
responsible for the evolution of the PDFs and do not contribute to the coefficient functions. 
In fact, in the variable-flavour- number scheme, they are MS-subtr acted from the coefficient 
functions and absorbed into the evolution of the PDFs. Thereafter, the asymptotic large- 

dependences of Ri at NLO should be proportional to ^^(Q^) and thus decreasing. This 
is familiar from the Callan-Gross ratio R = Fl/{F2 — Fl)^ as may be seen from its (x,(5^) 
parameterizations in Ref. [35j. Fortunately, this large-Q^ problem does not affect our results in 
Table [1] because the bulk of the HI data is located in the range of moderate values. 

The ratio Rc was previously studied in the framework of the /ct-factorization approach [29] 
and found to weakly depend on the choice of unintegrated gluon PDF and to be approximately 
X independent in the low-x regime (see Fig. 8 in Ref. [29 ). Both features are inherent in our 
approach, as may be seen at one glance from Eq. ([T6|) . The prediction for Rc from Ref. [29], 
which is included in Fig. [2] for comparison, agrees well with our results in the lower range. 
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10° 10^ 102 10^ 10^ 10^ 10^ 

GeV^ 

Figure 2: Rc, Rb, and Rt evaluated as functions of at LO from Eq. ([2T]) (dot-dashed lines) 
and at NLO from Eq. (p!6|) with /i^ = 4mf (dashed lines) and /i^ = + 4mf (solid lines). For 
comparison, the prediction for Rc in the /ct-factorization approach (dot-dot-dashed line) [29] is 
also shown. 

which supports the notion that the /c^-factorization approach partially accounts for the higher- 
order contributions in the low-x regime. 

7 Conclusions 

In this paper, we observed a compact formula [T3] for the ratio Ri = FI/F2 of the heavy- flavour 
contributions to the proton structure functions F2 and Fl valid through NLO at small values 
of Bjorken's x variable. We demonstrated the usefulness of this formula by extracting F2 and 
F2 from the doubly differential cross section of DIS recently measured by the HI Collaboration 
[21 [3] at HERA. These results agree with those extracted in Refs. ^ |3j well within errors. In 
the range probed by the HI data, NLO predictions agree very well with the LO ones and 
are rather stable under scale variations. Since we worked in the fixed- flavour- number scheme, 
our results are bound to break down for 4mf , which manifests itself by appreciable QCD 

correction factors and scale dependences. As is well known, this problem is conveniently solved 
by adopting the variable-flavour-number scheme, which we leave for future work. Our approach 
also simply explains the feeble dependence of Ri on x and the details of the PDFs in the \ow-x 
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regime. 
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